Black Hole Spectroscopy: Determining Waveforms from 3D Excited Black Holes 
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We present the first results for Cauchy nonlinear evolution of 3D, nonaxisymmetric distorted black 
holes. We focus on the extraction and verification of 3D waveforms determined by numerical relativ- 
ity. We show that the black hole evolution can be accurately followed through the ringdown period, 
and comparing with a recently developed perturbative evolution technique, we show that many 
waveforms in the black hole spectrum of modes, such as £ = 2 and I = 4, including weakly excited 
nonaxisymmetric modes (m 7^ 0), can be accurately evolved and extracted from the full nonlinear 
numerical evolution. We also identify new physics contained in higher modes, due to nonlinear 
effects. The implications for simulations related to gravitational wave astronomy are discussed. 
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Numerical relativity is the main tool for computing 
waveforms needed to aid in the detection and interpre- 
tation of gravitational wave signals from merging black 
holes (BH's). Recent and very detailed analysis shows 
that BH mergers may well be the first events seen, but 
that information from full 3D numerical simulations will 
be crucial in enhancing the rate of detection and in ex- 
tracting astrophysical parameters from the signals |jj . 

In response to this urgent need for 3D simulations, 
many groups worldwide have organized an effort to com- 
pute waveforms for binary BH mergers. Progress has 
been slow due to the difficulties of evolving the full Ein- 
stein equations, particularly in the BH case. For a review 
of this difficult problem, see Notable recent progress 
includes 3D characteristic evolution, success in apparent 
horizon boundary conditions (AHBC) and extrac- 

tion of waves from 3D simulations of axisymmetric BH's 
||. But, to date there have not been any published at- 
tempts to evolve, and extract meaningful waveforms from 
true 3D distorted BH's, as will be formed from the inspi- 
ral and merger of two colliding BH's, using full nonlinear 
numerical relativity. This ability is crucial for gravita- 
tional wave astronomy, and in this Letter we show that 
it can now be performed with a high degree of accuracy. 

We consider here the first detailed, nonlinear evolu- 
tions of nonaxisymmetric BH's, using the full machinery 
of standard 3+1 numerical relativity. We focus on a care- 
ful analysis of the 3D waveforms extracted, with an inde- 
pendent check on their accuracy and reliability, building 
particularly on two recent papers. Ref. || showed that 
one can use a full 3D nonlinear evolution code to study 
the excitation and ringdown of axisymmetric distorted 
BH's, including the extraction of the waveforms, such as 
those created by the head-on collision of two BH's. It was 
further shown, by comparison to 2D codes, that even very 
low amplitude waves could be extracted accurately in 3D. 



Ref. II developed a technique to use perturbation theory 
to evolve distorted BH data sets as a testbed for numer- 
ical relativity, showing how well this technique can work 
on a series of simulations performed with an axisymmet- 
ric distorted BH (2D) code, which allows one to verify 
the fully nonlinear code, and to identify nonperturbative 
effects. 

In this paper, for the first time we combine full 3D nu- 
merical relativity with perturbative techniques to show 
the extent to which nonaxisymmetric, distorted BH's can 
be evolved, and waveforms extracted accurately. We find 
that even very low energy [E < 10 _6 M, where M is 
the BH mass) , nonaxisymmetric wave modes can be reli- 
ably extracted in full 3D simulations. As emphasized by 
JlJ, nonaxisymmetric modes, particularly the £ = m = 2 
mode, may be the dominant signal in the final merger 
waveform. We show that not only can this mode be ac- 
curately extracted (albeit for nonrotating spacetimes at 
present), but even the much weaker I = 4, m family can 
be accurately computed. For a BH formed during some 
violent process (such as collision with another BH), we 
expect many modes to be excited. The main points of 
this paper are that (a) entire families of waveforms from 
the spectrum of an excited BH can now be reliably de- 
termined in full 3D numerical relativity, verified by per- 
turbation theory, and (b) nonlinear physics can also be 
identified. The combination of 3D nonlinear evolution 
and perturbative treatment is very powerful, and pro- 
vides an important testing ground for all 3D codes used 
to study waveforms from BH merger events, that will be 
crucial for gravitational wave astronomy. 

3D Simulations and Comparison with Pertur- 
bative Evolution. We consider the nonlinear numerical 
evolution of full 3D distorted BH's with a code written 
in 3D Cartesian coordinates, described in Refs. P,[7[-^0|. 
The code is built on earlier 3D work, where sphericaljg] 
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and colliding BH's Hi, and pure waves [0 were stud- 
ied, and was used in [[5||7]]^] , where axisymmetric distorted 
BH data sets were studied in 3D. Preliminary attempts 
to extract waveforms from 3D BH's were made in jf). 
This code solves the full set of Einstein equations with 
general gauge and slicing conditions, but here we use zero 
shift and the "1+log" algebraic lapse condition ||. 

The inner boundary on the BH is provided by an isom- 
etry condition H0,|l. Ultimately, one expects to replace 
this with an AHBC [pf, preventing the need to bend time 
slices so much that they ultimately destroy the simula- 
tions ||. The outer boundary conditions keep all func- 
tions fixed in time. If the boundary is far enough away 
this condition is usually adequate, but at present in 3D, 
without adaptive meshes (AMR) Jl2|,|l3|], it is generally 
impractical. This can create reflections of the waves 
at late times. Better boundary conditions, such as the 
perturbative outer boundary [0 , Cauchy-Characteristic 
matching [fl5|| , other treatments provided by hyperbolic 
formulations of Einstein's equations, or AMR should im- 
prove this situation. But these issues do not affect the 
main point we make below, which is that high resolu- 
tion and verifiable studies of truly 3D BH's, including the 
ringdown and propagation of very low amplitude nonax- 
isymmetric radiation modes, are now possible. 

The initial data we evolve in this paper belong to the 
family of distorted BH's discussed in Refs. I|]j]||[l6| pi , 
and correspond to a "Brill wave" jl9| superimposed on a 
BH. Such data sets mimic the endstate of two BH's collid- 
ing, forming a useful model for studying the late stages of 
BH mergers. To summarize we simply write the 3-metric 
here as d£ 2 = ip 4 (e 2q (drj 2 + d8 2 ) + sin 2 6W0 2 ) , where n 
is a radial coordinate related to the Cartesian coor- 
dinates by \J x 2 + y 2 + z 2 = e 1 ' ||. Given a choice 
for the "Brill wave" function q, the Hamiltonian con- 
straint leads to an elliptic equation for the conformal 
factor ip. The function q represents the wave sur- 
rounding the BH, and is chosen to be q (77, 9, (j>) = 

asm 11 9 fe~(^) +e~(^) ) (l + c cos 2 if). If the am- 
plitude a vanishes, the undistorted Schwarzschild solu- 
tion results; small values of a correspond to a perturbed 
BH. Although these data sets can include angular mo- 
mentum |^of , we focus here on the time symmetric case, 
so only the Hamiltonian constraint need be solved for 
the conformal factor i/j. Note that the form of the 3- 
metric, although fully nonaxisymmetric, does have a dis- 
crete (quadrant) symmetry about the z— axis and the 
equatorial plane. Hence we can evolve the BH in an oc- 
tant and save on the memory and computation required. 

Full 3D Nonlinear Evolution. The initial data set 
(a = —0.1,6 = 0, c = 0.5, w = l,n = 4) was evolved 
with the 3D code described above with a grid spacing 
of Ax = Ay = Az = 0.093M, with 300 grid zones in 
each coordinate direction. As in previous spherical and 
axisymmetric studies, the evolutions can be carried out 
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FIG. 1. We show waveforms for the I = m = 2 nonax- 
isymmetric mode extracted from the full nonlinear 3D simula- 
tion of the distorted BH. The top graph shows the waveform 
extracted at different locations, while the bottom shows the 
phase-shifted waveforms, allowing comparison between them. 

through about t — 40M, when the large gradients in met- 
ric functions, created by the use of singularity avoiding 
time slicings, cause the code to become rather inaccu- 
rate and crash. The problems encountered, not relevant 
for the main results of this work, should be eliminated 
through the use of AHBC (see, e.g. pi and references 
therein) . 

During the evolution, we extract the waves using a 
gauge-invariant technique developed originally by Abra- 
hams J2l],|2^] . The details of the present 3D application 
are given in [|)-||. Essentially, the Zerilli function ip is 
computed at various radii away from the distorted BH. In 
this work we extracted waveforms at Schwarzschild radii, 
r, of 8.0M, 10.2M, 12.6M and 14.9M. The 3D code and 
waveform extraction procedure were developed and veri- 
fied through careful comparison to an axisymmetric code 
|l6| , ^7| and also on full 3D initial data || . Here we apply 
the same techniques to study the first waveforms emitted 
from nonaxisymmetric distorted BH's in full numerical 
relativity. 

In Fig. [l] we show the nonaxisymmetric i — m = 2 
mode emitted during the evolution of this distorted BH. 
Waves are extracted at the different radii, showing the de- 
velopment of this mode in time. The quasinormal ringing 
of the BH is clearly evident and has the correct complex 
frequency. The modes are normalized so that the energy 
they carry is given by E — (l/327r)-0 2 . For this mode, 
we find that the energy is about 10 _6 M. 

As mentioned, this 3D mode is expected to be an im- 
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FIG. 2. We show the waveform for the £ — 2, m = mode, 
extracted from the linear and nonlinear evolution codes. The 
dotted (solid) line shows the linear (nonlinear) evolution. 

portant and dominant mode in the spectrum of BH ex- 
citations created during the merger process. This is the 
first example showing that nonaxisymmetric modes can 
be computed in a full 3D simulation, in spite of all the dif- 
ficulties of BH simulations. The calculation was repeated 
at higher and lower grid resolutions, the waveform was 
found to be converging to between first and second order 
in the grid spacing, with the changes in the waveform 
at different resolutions being minimal. Beyond resolu- 
tion studies, having an independent technique to deter- 
mine the expected waveform, to which comparisons can 
be made, is also essential in studying their quality. As 
emphasized elsewhere perturbative techniques can 
also be used as a testbed to study black hole dynamics, 
and provide physical insight into the nature of the solu- 
tion. We turn to perturbation theory in the next section 
to provide this testbed, and the physical insight it brings. 

Comparison with Perturbative Evolution. Although we 
have seen that the 3D modes have the right quasinor- 
mal frequency, without a careful check of the simula- 
tions, we have no guarantee that the results are actually 
correct. However, building on prior work of Abrahams, 
Price, Pullin and others (see |^3) for recent review), we 
have recently developed a technique to study distorted 
BH's through a perturbative approach to the evolution 
H . The technique makes use of the extraction technique 
used above to extract initial data for the Zerilli equation, 
describing linearized perturbations of spherical BH's. As 
long as the initial wave content on the Schwarzschild 
background is small, one should be able to use this sim- 
ple linear wave equation to evolve the gravitational waves 
and predict the waveforms that must be seen by the full 
3D, nonlinear code. We have taken such initial data, 
evolved it with the Zerilli equation, and compared with 
waveforms extracted from the full 3D nonlinear evolution 
in Cartesian coordinates. 

In Figs. |[ H, and || we show results for a large part of 
the BH spectrum of modes excited. The waveforms for 
the linear and nonlinear evolutions are each plotted on 
the same graphs, extracted at r = 12. 6M. The total en- 
ergy radiated for these modes runs from £~3x 10~ 4 M 
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FIG. 3. Waveforms are shown for the I = 4, m — 
extracted from the linear and nonlinear evolution codes. The 
dotted (solid) line shows the linear (nonlinear) evolution. 

for the £ = 2, m = mode, to £ ~ 3 x 1CT 7 M for the 
I = 4, to — 2 mode. The agreement between these two 
completely independent treatments is remarkable, giving 
complete confidence in the reliability of these waveforms. 
Not only do the perturbative results confirm those of full 
3D numerical relativity, but the 3D results confirm the 
perturbative treatment. This is an important point to 
emphasize, as it shows that these modes are actually in 
the perturbative regime. Hence the gravitational wave 
physics of these modes is adequately treated by linearized 
theory. This must be captured by the full 3D simulation, 
even though the singularity avoiding slicing chosen for 
the 3D simulation forces the BH background to evolve 
in a fully nonlinear manner (leading to "grid stretching" , 
etc.) The £ — to = 2 waveform already shown above has 
the same level of agreement with its perturbative coun- 
terpart. A full analysis of these waveforms and other 3D 
BH spacetimes will be published elsewhere Q. 

We emphasize that although the amplitude of distor- 
tion of this BH is low enough that some modes can be 
treated with linear perturbation theory, there are lim- 
its to its applicability. Other modes can be out of the 
linear regime, and nonlinear evolution should not agree 
with linearized evolution results. Detailed analysis of the 
initial data shows that only the modes I = 2, to = 0, 2 
and £ — 4, m = 0, 2 should occur at linear order in 
the wave amplitude a, all other modes occurring at or- 
der a 2 or above. These higher order modes cannot be 
accurately evolved with the first order Zerilli equation; 
Ref. [£5| shows that in such cases there is a source term 
arising from nonlinear mixing of the modes that appear 
at linear order. In some sense, these modes are too small 
for linear treatment to adequately capture the physics, 
yet as we show now they can be treated in full 3D nu- 
merical relativity. 

In Fig. ||, we consider such a low amplitude case, I — 6, 
to = 2. (Other modes are more difficult to extract be- 
cause of numerical resolution, e.g. £ — to = 4 has not 
converged in our simulations. Full details will be pub- 
lished elsewhere). This is a difficult mode to extract nu- 
merically in 3D p4] , but the waveform is not significantly 



3 



0.002 - 



-0.003 - 




10 20 30 40 

Schwarzschild time, t/M 

FIG. 4. Waveforms are shown for the I = 4, m = 2 mode, 
extracted from the linear and nonlinear evolution codes. The 
dotted (solid) line shows the linear (nonlinear) evolution. 



in gravitational wave astronomy. 
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FIG. 5. Waveforms are shown for the £ = 6, m = 2 mode, 
extracted from the linear and nonlinear evolution codes. The 
dotted (solid) line shows the linear (nonlinear) evolution. The 
discrepancy is attributed to a nonlinear effect. 



altered by more resolution, and show a marked disagree- 
ment in contrast to the previous Figs. The linear energy 
calculation indicates that the energy of this mode is only 
E ~ 10~ 10 M . Simulations at different values of a show 
the amplitude of this waveform scaling as a 2 . Hence, 
we fully expect the perturbative approach to fail for this 
mode, and we must turn to nonlinear treatments, such as 
a 3D numerical relativity code, to compute it. This same 
effect is discussed in || for highly resolved 2D cases where 
nonlinear mode mixing is very clearly seen; details of the 
more difficult 3D cases are in preparation for publication 
elsewhere p4| . 

Conclusions. We have shown for the first time that 
evolutions of a new class of full 3D distorted black 
holes can be accurately performed in 3D Cauchy non- 
linear evolution codes in Cartesian coordinates, and fur- 
ther that a large spectrum of different nonaxisymmetric 
waveforms can now be determined from the simulations, 
even when they carry a tiny fraction of the ADM mass 
(E ~ 10 _6 M). Such capability will be crucial in simulat- 
ing 3D black hole mergers, needed for effective detection 
of waves. We have also applied newly developed pertur- 
bative techniques to verify that the waveforms are very 
accurate, and also to show nonlinear effects. This com- 
parison technique should provide a powerful testbed and 
analysis tool for all 3D codes under development for use 
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